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Recent measurements of the London penetration depth tensor in the cuprates find a weak tem- 
perature dependence along the c-direction which is seemingly inconsistent with evidence for d-wave 
pairing deduced from in-plane measurements. We demonstrate in this paper that these disparate 
results are not in contradiction, but can be explained within a theory based on incoherent quasipar- 
ticle hopping between the Cu02 layers. By relating the calculated temperature dependence of the 
penetration depth Ac(T) to the c-axis resistivity, we show how the measured ratio \c(0) / Xc{T) can 
provide insight into the behavior of c-axis transport below Tc and the related issue of "confinement." 



PACS numbers; 74.20.-z. 74.25.-q, 74.25.Nf, 74.25.Fy 



Measurements of the temperature dependence of the 
in-plane penetration depth in YBa2Cu307_5 (YBCO) 
have been interpreted as-.providing strong support for a 
d-wave order parameter .B With the availabihty of high- 
quahty single crystals, these measurements aijcjiow be- 
ing extended to the c-axis direction in YBCOHa as well 
as in BjUiSrsCaCusOs (BSCCO)!^ and Lai_xSrxCu04 
(LSC0).l3 There appears to be one consistent feature of 
these different experiments: the (low) temperature de- 
pendence of the c-axis penetration depth Ac(T) is much 
weaker than that observed in the ab plane. In addition, 
the penetration depth ratio A^(0)/A^(r) in YBCO is lin- 
ear in T at low T with a small slope that decreases with 
decreasing oxygen content 

The present paper addresses these penetration depth 
data in conjunction with the c-axis resistivity Pc{T) for 
a variety of different cuprates. We argue that the tem- 
perature dependence of Pc{T) essentially determines that 
of Xc{T)\ in particular, \c{T) is expected to have a weak 
temperature dependence at low T whenever Pc{T) ex- 
hibits a strong semiconducting behavior. In this way, we 
can reconcile the a6-plane data supporting d-wave pair- 
ing with the observed c-axis behavior. In addition, by 
presenting typical results for Ac with the corresponding 
Pc we illustrate the generic features of these quantities 
without relying on specific parameterizations of (or fits 
to) existing data. Through the study of the c-axis cou- 
pling, we touch on the issues which may lie at the heart 
of the nature of the normal stateQ and also potentially 
the mechanism of high-temperature superconductivity.Q 
In particular, our results imply that a Fermi-liquid-based 
description of the electronic states in the Cu02 planes 
appears as consistent with the data as theories based on 
the idea of "confinement," even though the underlying 
assumptions are considerably different. 

To understand the behavior of \c{T) and the nature of 
c-axis coupling below , we utilize an incoherent hopping 
model for the c-axis couplingJjO Above the critical tem- 
perature Tc, this model views each Cu02 layer as a two- 



dimensional Fermi liquid which is weakly coupled to its 
nearest neighbors. This model yields a semiconducting 
T dependence of Pc{T) idpciT) / dT < 0) while maintain- 
ing a finite residual resistivity Pc(0). The semiconducting 
behavior is also associated with the absence of a Drude 
peak in the c-axis optical conductivityH Thus, while the 
c-axis transport properties suggest an insulating state at 
low temperatures, the actual zero-temperature state is 
nevertheless metallic.B Extending this model to Ac and 
comparing its predictions with the available experimen- 
tal data provides an opportunity to learn about the pre- 
cise low-T behavior of the c-axis resistivity as well as the 
degree to which the layers communicate in the supercon- 
ducting phase. 

We begin by writing the Hamiltonian for the electronic 
system as the sum of a Hamiltonian for the individual 
Cu02 layers, X^m^mi inter layer coupling term 

H±: Hc\ = J2m^"^ + We leave the intra-layer 

Hamiltonian unspecified except to demand that each Hm 
yield a two-dimensional Fermi liquid which becomes a ei- 
ther an s- or a d-wave superconductor below a critical 
temperature Tc- H± is taken to beETLJ 



(1) 



where = + Vi„i + j 9i—3-m4'jra and Cima is the 
usual quasiparticle annihilation operator for site i in 
layer m and spin projection a. Physically, the interlayer 
coupling arises from quasiparticle hopping due to wave- 
function overlap (parameterized by tj^), impurity scatter- 
ing (modeled by the random variable 1^™), and bosonic 
scattering (represented by the field 0™ which couples to 
the electronic quasiparticles with strength gim)- 

In treating _ff_L, we are guided by two observations: 
(1) the mean free path in the c-direction extracted 
from normal-state .tran-sport measurements is less than 
the lattice spacing ,ll3ii3 and (2) the c-axis properties in 
the superconducting state are consistent with a picture 
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where nftarest-.neighbor Cu02 layers form an SIS tunnel 



junctionD'liilo These observations suggest that the c-axis 
transport may be viewed as an incoherent tunneliga 
cess. Several theories of this incoherence exist,LfHll2 
but for simplicity we adopt the phcnomenological model 
of Rcf. H and simulate the effects of incoherence by per- 
forming our calculations to second order in H± . This prpn 
cedure immediately yields an intrinsic Josephson effecttHI 
and a reasonable piaEnitude and T dependence of the 
c-axis resistivityj3~lij'Lj 

In the normal state, we can calculate the c-axis dc con- 
ductivity by direct analogy with the problem of tunnelina 

3nim 



in a normal metal-insulator-normal metal junctior 
and obtain a c-axis conductivity which is the sum of the 
conductivities due to each process: CTc = o-d"°ct _|_ ^imp ^ 
^mci|^ Generalizing the results of Ref. |l^ to anisotropic 
scattering, we can write these terms as 



^direct 
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Here, A^o is the density of states per unit cell per spin at 
the Fermi surface, Ti/rab is the intra-layer scattering rate, 
(To = (47re^/?i)(o?/a^), d (a) is the inter- (intra-) layer unit 
cell dimension, e is the electronic charge, and the an- 
gle brackets denote a normalized average over the Fermi 
surface. For simplicity, we take the interlayer inelastic 
scattering to have an Einstein spectrum with frequency 
fig but keep the wave- vector dependence in the impurity- 
(V^) and boson-assisted [g^) hopping amplitudes. 

The calculation of the penetration depth in the su- 
perconducting state may be performed either by not- 
ing that the c-ajds critical current cx in layered 
superconductorsEiapd computing the critical current in 
the standard way,E3'E3 or by calculating the optical con- 
ductivity and extracting the penetration depth from the 
imaginary part.cil In either case, one obtains 



where c is the speed of light, T is the temperature (fcs 
= 1), iV|| is the number of lattice sites in a single layer. 
ki = (k, iw;), F{ki) is the Gor'kov propagator ,E2lEJ and 
t^{ki, k[i) is the generalized interlayer hopping matrix el- 
ement. As with the conductivity, this matrix element 
is the sum of the matrix elements due to direct scatter- 
ing, fc;',) — ^(5k,k"5i,;' t\; impurity-assisted 
scattering, t'^'^™P{ki, k'l,) — ^6ij' |Vk_k'|^; and boson- 
assisted scattering, P'^'^°^{ki, kl,) = |5k-k' |^ -D(fc; — k'l,). 



In the last term, D{ki) is the propagator for the inter- 
layer boson. We note that all three mechanisms pro- 
duce a contribution to the penetration depth. This re- 
sult is not surprising, given the close similarity between 
incoherently coupled layers and tunnel junctions: the 
direct and impurity-assisted hopping processes formally 
resemble lliQ— iwcesses considered in conventional SIS 



junctions J 



and boson-assisted hopping has been 



known to contribute to tunneling iii superconductim 
junctions in both the quasiparticleCJ and JosephsorO 
channels for some time. 

To compute the penetration depth from these formu- 
lae, we make several simplifying assumptions. First, we 
use the standard BCS form for the Gor'kov propagators 
and perform the sums over k in the usual way by re- 
stricting the wave vectors of the self-energies and matrix 
elements to the Fermi surface and then integrating the re- 
maining energy dependence from —oo to +oo. We reiter- 
ate that the weak interlayer coupling allows us to perform 
our calculations to second order in the interlayer hop- 
ping amplitudes, and this implies that inter-layer scatter- 
ing effects which act ia reduce 1/A^ are of higher order 
and can be neglected. Second, we make the reasonable 
approximation that the layers are identical. Third, we 
account for the anisotropy of the hopping processes b;^ 
taking the impurity scattering matrix element |Vk_k'| 
to have a Lorentzian form with a maximum at k = k' 
and a half width 5k/kp — 0.01 and by assuming that the 
wave- vector structure of Igk-k' | is such that the boson- 
assisted processes contribute to both the resistivity and 
the penetration depth with the same strength. The qual- 
itative trends in pc and Ac we discuss do not depend on 
these choices; see Ref. |ll| for details. 

From these relations, we see that pc — I/cfc is de- 
termined by the same parameters as Ac. We are there- 
fore able to connect the two quantities and examine the 
qualitative predictions of our theory. Fig. |l| shows the 
resulting curves for Pc(T) and the corresponding c-axis 
penetration depth ratio for both s- and d-wave pairing in 
several limiting cases of our model. 

If only direct hopping is present [cf. Fig. |l|(a)], only 
those terms corresponding to wave function overlap con- 
tribute to the resistivity and the penetration depth. This 
case may be relevant for materials like fully oxygenated 
YBCO, which is one of the least anisotropic cuprates 
{pc{Tc)/pab{Tc\^ 10-102 compared to Pc{Tc) / Pab{Tc) ~ 
10^ for BSCCOIIj) and therefore potentially the least in- 
coherent. In this instance, pc reflects the temperature 
dependence, though not the magnitude, of the a6-plane 
resistivity: pc oc pab « T [cf. Eq. (||)]. The resulting 
penetration depth is also marked by this near-coherence 
and has the same temperature dependence for both the 
ab- and c-axis directions, regardless of pairing symmetry. 

The case of assisted hopping, where only the impurity- 
and boson-assisted processes contribute to the c-axis 
transport, is shown in Fig. |l|(b). Materials like BSCCO 
with c-axis mean free paths much less than the lattice 
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spacing are expected to be close to this limit, pc in 
this case is marked by a negative temperature derivative 
which gives rise to a semiconducting temperature depen- 
dence above T^. This upturn in the resistivity is due to 
the freezing-out of the inelastic interlayer scattering at 
low T, which inhibits c-axis transport and therefore in- 
creases the resistivity. At lower T, however, the impurity 
scattering acts to limit the conductivity, and Pc{T) sat- 
urates. The c-axis penetration depth is also modified by 
the incoherent transport and becomes distinct from the 
a6-plane result. For both pairing symmetries, the c-axis 
penetration depth ratio is larger than the ah plane pen- 
etration depth ratio at all temperatures. In particular, 
the d-wave penetration depth ratio A^(0)/A^(T) resem- 
bles the s-wave case, although the low-temperature be- 
havior is still a power law. In addition, the difference 
between in-plane and c-axis penetration depth ratios is 
much more pronounced for d-wave pairing than for s- 
wave paring (see the upper inset in Fig. 0(b)). Moreover, 
the temperature dependence in the d-wave case is close 
to, and may be experimentally indistinguishable from, 
the (s-wave) Ambegaokar-BaratofS form. 

In the intermediate case, all processes contribute to 
c-axis transport [Fig. 0(c)]. Compounds such as de- 
oxygenated YBCO may belong to this class. As in the 
assisted hopping case, pc looks semiconducting at high 
temperatures. At low temperatures, however, the in- 
plane scattering rate is reduced, and this leads to a low-T 
conductivity dominated by direct hopping processes [cf. 
Eq. (H)]. The net result is a peak in p^ which may lie 
below Tc- In contrast to the resistivity, the penetration 
depth ratio shows no new behavior in this case, but is 
mid-way between the direct result, where both ah- and 
c-axis penetration depth ratios are the same, and the 
assisted hopping result, where they are considerably dif- 
ferent. Note especially that the low-temperature c-axis 
penetration depth ratio for d-wave pairing is clearly lin- 
ear in T but with a much smaller slope than its in-plane 
counterpart. 

Fundamentally, the behavior discussed above results 
from the different characteristic temperature depen- 
dences that arise when the interlayer coupling is very 
weak (incoherent) as opposed to when it is strong (co- 
herent). In a material with very weak interlayer cou- 
pling, quasiparticle transport between adjacent layers is 
analogous to tunneling in NIN or SIS junctions. This 
tunneling is mediated by both elastic are inelastic scat- 
tering and is therefore associated with a semiconduct- 
ing c-axis resistivity in the normal state and a very low 
slope of the c-axis penetration depth ratio at low tem- 
peratures. The small low-tcmperature slope of X^'^{T) 
is a special feature |Of specular Josephson tunneling and 
occurs in either s— Ej or d-wavdiil superconductors. On 
the other hand, in materials with stronger interlayer cou- 
pling, quasiparticle transport is nearly coherent and so 
quasi-three-dimensional results obtain: pc is metallic, 
and A^(0)/A^(T) possesses a larger low-T slope for d- 
wave pairing relative to the weakly coupled case. 



Empirically, therefore, one should associate a reduc- 
tion in the low-T slope of the c-axis penetration depth 
ratio in either s- or d-wave superconductors with a de- 
crease in the interlayer coupling and a semiconducting 
temperature dependence of pc- Our theory thus accom- 
modates the currently available data in both the intra- 
and interlayer directions despite their very different low- 
temperature slopes. Moreover, this consistency suggests 
that the cuprates may indeed be d-wave superconductors. 
While our theory can explain the qualitative features of 
the present experimental data, these data are incomplete; 
further systematic experiments on different cuprates and 
for particular cuprates with different stoichiometrics are 
clearly called for to further test the trends reported in 
this paper. In addition, further theoretical effort is re- 
quired to understand the relationship of our theory to 
others in the literature and to provide a way of distin- 
guishing them experimentally. 

Among these other tlaeories of c-axis coupling, the 
"confinement" approachQ is worth discussing further, 
since we come to similar conclusions despite vastly differ- 
ent starting assumptions. In contrast to our Fermi-liquid- 
based approach, the confinement theory asserts that each 
Cu02 layer is, a spin-charge separated tomographic Lut- 
tinger liquid.El Nevertheless, the expressions for the inter- 
layer transport which arise in this model are similar to 
those in this work with two important differences: (1) the 
Green's functions used in the calculation correspond to a 
Luttinger and not a Fermi liquid, and (2) inelastic scat- 
tering is an intrinsic part of the interlayer transport pro- 
cess and is not extrinsic to the layers as in, for example, 
the boson-assisted hopping term in our theory. Despite 
these differences, bath theories obtain a semiconducting 
T dependence of p^ aad a small low-T slope of the pen- 
etration depth ratiOjCB in agreement with experiments. 
This comparison suggests that models of c-axis coupling 
in the cuprates should view interlayer transport as a tun- 
neling process with a strong inelastic component, but the 
detailed origin of this effect may not be readily extracted 
from the currently available experiments. 
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FIG. 1. Penetration depth ratios Xlt{0) / >^lb{T) (solid 
line) and X^{0)/\^{T) (dot-dashed line) as a function of the 
normalized temperature T/Tc for both d-wave (main figure) 
and s-wave (upper right inset) pairing in different limits of the 
incoherent hopping model described in the text. The corre- 
sponding c-axis resistivity pc is shown in the lower-left inset 
normalized to its value at Tc as a function of temperature. 
The limits illustrated are (a) direct hopping, which may be 
relevant for YBa2Cu306.9, (b) assisted hopping, which may 
be the case in Bi2Sr2CaCu2 08, and (c) a combination of these 
processes, which may obtain in underdoped compounds like 
YBa2Cu306.4. Note in (a) that both the ab- and c-axis pen- 
etration depth ratios are the same when only direct hopping 
is present. 
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